ASPECTS OF TOEPLITZ DETERMINANTS 

I. KRASOVSKY 

Abstract. We review the asymptotic behavior of a class of Toeplitz (as well as related 
Hankel and Toeplitz + Hankel) determinants which arise in integrable models and other 
contexts. We discuss Szego, Fisher-Hartwig asymptotics, and how a transition between them 
is related to the Painleve V equation. Certain Toeplitz and Hankel determinants reduce, in 
certain double-scaling limits, to Fredholm determinants which appear in the theory of group 
representations, in random matrices, random permutations and partitions. The connection 
to Toeplitz determinants helps to evaluate the asymptotics of related Fredholm determinants 
00 ■ in situations of interest, and we review the corresponding results. 

O-t- 1. Introduction 

Let f(z) be a function integrable over the unit circle C with Fourier coefficients 

fj= ^r f( e ^ e ~ ijede > * = o, ±i, ±2, . . . 



O 



Then the n- dimensional Toeplitz determinant of a Toeplitz matrix with symbol f(z) is given 



00 ! by 

<N ' (1-1) D n (f) = detC/;,-*)^. 

Substituting here the expressions for the Fourier coefficients, and using formulae for Vander- 
q ■ monde determinants, one obtains another useful representation: 

•I nl-K r-2-E n 

(!- 2 ) D ^ = TM ■■■/ II \e^-e^Uf(e^)d6 r 

V ' ' J0 J0 l<j<k<n j=l 

Toeplitz determinants are closely related to the polynomials orthogonal with weight f(z) 
on the unit circle. Namely, if -Dfc(/) ^ 0, k = k Q , k + 1 . . . , for some ko > 0, then the 

polynomials <f>k{z) = Xk% k + ■ ■ ■ , 4>k( z ) = Xk zk H of degree k, k = k ,k + l, . . . , satisfying 

(1-3) ^j* <p k (z)z-if{z)d6 = x" 1 ^, -LJ^^z-iyf^de = x^Sjk, 

z = e i0 J j = 0,l,...,k, 

exist and Xk = yDk/Dk+i, where, by convention, D = 1. We see from (1.2) that if f(z) is 
positive on C, we have D n (f) > for all n, and therefore in this case we can set k = 0. 

A Toeplitz determinant can be represented as a Fredholm determinant of an integral 
operator acting on L 2 (C) which belongs to the special class of so-called integrable operators 
[48]. It can also be written in a different way in terms of a Fredholm determinant of an 
operator now acting on £ 2 (n, n + 1, . . . ) [62, 27, 30]. (Note that the symbols f(z) considered 
in [62, 27, 30] are assumed to be sufficiently smooth.) Another useful property of many 

l 
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.D n (/)'s encountered in applications is the existence of simple differential identities relating 
the determinant to orthogonal polynomials evaluated at a few special points. The precise 
form of such identities depends on the given f(z). 

It turns out that the above properties play a key role in making Toeplitz determinants 
amenable to a detailed asymptotic analysis, in particular, by Riemann-Hilbert-Problem 
methods. 

In this paper, we will review some asymptotic results on D n (f) (and related Hankel, 
Teoplitz+Hankel, and Fredholm determinants) and briefly mention their applications in 
integrable models, random matrices, random permutations, group representation theory, 
and also in various conjectures on Riemann's ( and Dirichlet's L-functions. This review is 
based, to a large extent, on the recent work of the author with T. Claeys, P. Deift, A. Its, and 
J. Vasilevska. For other aspects of Toeplitz determinants not mentioned here, the reader is 
referred to [32, 34, 33, 67, 47, 49, 50] for properties of Toeplitz matrices and determinants, to 
[34, 108, 109, 110, 89, 64, 98, 19, 91, 90, 105, 72, 69, 14] for generalizations to the continuous, 
higher-dimensional, and the block- Toeplitz cases (with relations to stationary determinantal 
processes, integrable models, and entanglement entropy), and to [88, 1, 20] for connections 
with multiple orthogonal polynomials. 

The paper consists of three parts: in Section 2, the simplest asymptotics with f(z) fixed 
and n — > oo are considered; in Section 3, the symbol f(z) is allowed to depend on n in 
ways which describe a transition between different asymptotic regimes arising in Section 2; 
in Section 4, the symbol f(z) also depends on n, but in such a way that in the limit n — > oo, 
Toeplitz determinants turn into certain Fredholm determinants which are important for 
random matrices and random permutations. Following standard practice, we refer to the 
large n asymptotics of Sections 3 and 4 as double-scaling limits. 

2. Asymptotics for a fixed symbol 

We assume in this section that f(z) does not depend on the size of the determinant n. 
We are interested in the asymptotics of D n (f) as n — > oo. 
The following result is basic. 

Theorem 2.1 (Strong Szego limit theorem). Let f(z) be non-zero on C, In f(z) G L 1 (C) 7 
and suppose that the sum 

(2.1) S(f)= J2 \ k W(^f)k\ 2 , (In f)k = T \nf(e ie )e- ike de, 

k=-oo n ^° 

converges. Then 

oo 

(2.2) In D n (t) = n(ln/) + ^ fc(ln/) fc (ln/)_ fc + o(l), as n ->■ oo. 

fc=i 

The theorem was initially proved by Szego [99, 100] (the leading term in 1915, the next 
in 1952) under stronger conditions on f(z). The conditions were then weakened by many 
authors. In the present form, the theorem was proved in [68, 65, 80]. See [97] for a detailed 
account. 

A strong motivation to study such asymptotics first came in the end of 1940's after On- 
sager's solution of the 2-dimensional Ising model and his observation that a 2-spin correlation 
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function in the model can be written as a Toeplitz determinant D n (f), where n denotes the 
distance between the spins. For temperatures less than critical (T < T c ), the symbol of this 
Toeplitz determinant has an analytic logarithm in a neighborhood of the unit circle and, 
moreover, (ln/) = 0. Therefore, Szego's theorem can be applied, and one concludes that 
D n (f) tends to a constant as n — > oo. Thus the correlation does not decay as the distance 
increases, which indicates the presence of a long-range order, and hence, a magnetization. 
As T /* T c , however, 2 singularities of f(z) approach the unit circle at z — 1, and, at T = T c 
merge into a single singularity on C; namely, a jump-type singularity at z = 1 (see, e.g., 
[94]). For f(z) with such a singularity, the sum (2.1) diverges, and Theorem 2.1 can no longer 
be applied. In fact, it turns out [94] that in this case D n (f) decays as rr 1 ^ and, therefore, 
there exists no long-range order. For correlation functions arising in other situations, such as, 
e.g., the so-called emptiness formation probability in the XY spin chain in a magnetic field 
[40, 61], one obtains Toeplitz determinants with both jump-type and root-type singularities, 
and in the most general situation one is led to consider symbols of the form: 

m 

(2.3) f(z) = e ^W=o& Y[ \z - z 3 \^g Zj ^(z)zJ^, z = e ie , 9 e [0, 2k), 

3=0 

for some m > 0, where 

(2.4) Zj = e*i, j = 0,...,m, = 6 < 6 1 < ■ ■ ■ < 6 m < 2ir; 



(2-5) g Z]A (z) 



& i < arg z < 9j 
e -iirpj q. < arg z < 2tc 



(2.6) 9taj > -1/2, pjEC, j = 0, . . . , m, 

and V(e ze ) is a sufficiently smooth function on the unit circle (see below) with Fourier 
coefficients 

(2.7) V k = ±- j" V{e i9 )e- ki9 d0. 

2vr Jo 

The canonical Wiener- Hopf factorization of e v ^ is given by 

(2.8) e v{z) = b + (z)e v %„(z), b+(z) = e^^ VkZ \ b_{z) = e^*=-« VkZ " . 

The condition (2.6) on aj ensures the integrability of /. Note that the size of the jump 
at Zj is determined by the parameter j3j, and the root-type singularity, by aj. We assume 
that Zj, j = 1, ... ,m, are genuine singular points, i.e., either aj ^ or f3j ^ 0. However, 
the absence of a singularity at z — 1, i.e. the case a$ = (3o = 0, is allowed. Singularities 
of type (2.3) are known as Fisher-Hartwig singularities because of the work [58] where the 
authors summarized a variety of applications of Toeplitz determinants with such symbols 
and presented a conjecture about the asymptotic form of D n (f) in this case. Due to the 
subsequent efforts of many workers, we have the following description of the asymptotics. 

Define the seminorm 

(2.9) IH0III =max|9fc&-&/9 fc |, 

where the indices j,k — are omitted if z — 1 is not a singular point, i.e. if «o — Po — 0. 
Note that in the case of a single singularity, we always have \\\/3\\\ = 0. 
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First, consider the situation when \\\/3\\\ is strictly less then 1. 

Theorem 2.2. Let f(z) be defined in (2.3), \\\(3\\\<l, ftctj > -1/2, a 3 - ± fa ^ -1, -2, . . . 
for j,k = 0,1,..., m, and V(z) satisfies the smoothness conditions (2.11), (2.12) below. 
Then as n — > oo ; 



(2.10) D n (f) =exp 



nV + J2 kV k V - 



fc=i 



j=0 

ctjP k -a k pj 



C\< -i<^b<rn \ 3 / 

" ^(1 + ^ + ^)^(1 + ^-^) 

11 r-n 4-9™/i i 1-1- ^ 1 ^ 



0<j'<fc<m 



m 
X 



G(l + 2 0j 



where G(x) is Barnes' G-function [10]. ITie double product over j < k is set to 1 i/ 
m = 0. T/ie branches in (2.10) are determined as follows: b±(zj)~ aj± ^ j = exp{(— a 3 ± 
fa) EZi V± k z ±k }, ( ZkZ -i e -™) a ^- a ^ = exp{i(6 k - 9j - n)(a 3 /3 k - a k fa)}. 



Note that since G(—k) =0, k — 0,1, ... , formula (2.10) no longer represents the leading 
asymptotics if a 3 + fa or a 3 — fa is a negative integer for some j. Such degenerate cases can 
be handled by carrying the analysis to higher order, but we present no further details here. 

The smoothness condition on V(z) assumed in Theorem 2.2 is that 



oo 



(2.11) 2^ \k\ s \V k \ < oo 

k=— oo 

holds for some s such that 

(2 12) s > i+_e.™ o k^,) 2 + tm 



i- 

Note that the condition \\\/3\\\ < 1 is important here. 

The Barnes' G-function first appeared in asymptotic Toeplitz theory in the work of Lenard 
[92]. Theorem 2.2 was proved by Widom [104] in the case when $ta 3 > —1/2, and all fa = 0, 
and with a stronger condition on V(z). In [11], Basor extended the result to tftotj > —1/2, 
$lfa = 0, and in [12], to a 3 = 0, \3lfa\ < 1/2. In [31], Bottcher and Silbermann established 
the result in the case that \$laij\ < 1/2, \$lfa\ < 1/2. In [53], Ehrhardt proved the theorem 
for the full range of parameters, namely ffiatj > —1/2, |||/3||| < 1, and for G°° functions 
V(z). These results were established by operator-theory methods (see [53] for a review of 
these and other related results including an extension to 9fta < —1/2, 2a ^ — 1, —2, . . . , 
when / is replaced by a suitable distribution). In [44], the authors reprove the theorem by 
Riemann-Hilbert-Problem methods, and relax the smoothness conditions on V(z) to (2.11), 
(2-12). 

Consider now the general case of Fisher-Hartwig symbols f(z) with the restriction \\\/3\\\ < 
1 removed. Note first that f(z) has several representations of type (2.3) with different sets 
of parameters fa. Namely, if each fa in (2.3) such that either fa ^ or a 3 ^ is replaced by 
fa = fa + rij, where n 3 are integers subject to the condition YT 3 =q n j = 0> then the resulting 
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function f(z; no, ... , n m ) is related to f(z) in the following way: 

f(z;n ,...,n m ) = J\zJ nj f(z), 

3=0 

i.e., it differs from f(z) only by a constant. Each f(z;n , . . . ,n m ) so obtained is called 
a FH-representation of the symbol. Denote by At the (finite) set of FH-representations 
for which Y^jLoffiftj) 2 * s minimal- There exists a simple procedure (see [43]) to solve this 
discrete variational problem and to construct At explicitly. One can show that there is 
always a FH-representation with |||/9||| < 1, and we have the following 2 mutually exclusive 
possibilities: 

• If there exists a FH-representation such that |||/3||| < 1 then it turns out that this 
FH-representation is the single element of At. In particular, if |||/3||| < 1, the set 
AA consists of a single element corresponding to all rij = 0, and Theorem 2.3 below 
reduces to Theorem 2.2. 

• If there exists a FH-representation such that |||/3||| = 1 then At consists of several 
(at least 2) elements. 

The set At is called non-degenerate if it contains no representations for which otj + f3j or 
aj — Pj is a negative integer for some j. The general result is as follows. 

Theorem 2.3. Let f(z) be given in (2.3), V(z) satisfy the condition (2.11) above for some 
sufficiently large s (depending only on aj, (3j), and 9ftaj > —1/2, f3j E C, j — 0, 1, . . . ,m. 
Let At be non- degenerate. Then, as n — > oo, 

(m \ n 

\{z n A K(f(z;n ,...,n m ))(l + o(l)) 

where the sum is over all FH-representations in At. Each lZ(f(z; no, . . . , n m )) stands for the 
right-hand side of the formula (2.10), without the error term, corresponding to f(z; n , . . . , n m ) 

An explicit lower bound on s (depending on /3j, aj) similar to (2.12) is given in [43]. 

This theorem was conjectured by Basor and Tracy [18] and proved in [43]. 

Hankel and Toeplitz+Hankel determinants are also of interest. Let w(x) be an integrable 
function on a subset J of KL Then the Hankel determinant with symbol w(x) supported on 
J is given by 

(2.14) D*(w(x)) = det ( / x J+k w(x)dx j 

\JJ / j,k=0 

When J is a finite interval - we then set J = [— 1, 1] without loss of generality - Hankel 
determinants are related to Toeplitz determinants by the following formulae [43], involving 
the orthogonal polynomials (1.3): 

f(e ie ) 

(2.15) W ( x ) = l±—L j x = cos6, se[-l,l]; 

| sin 0\ 

/ 01 m \riH/ I \\i2 n n \X2n + 02n(O)) / r / \\ 

(2.16) [D n (w(x))] = 4(w _ 1)a 02w(1)02ra( _ l} DMz)). 
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A particularly interesting class of Toeplitz+Hankel determinants appearing in the theory 
of classical groups and its applications to random matrices and statistical mechanics (see, 
e.g., [7, 60, 81]) is defined as follows for even /(e) = f{e~ ie ) (for even / the matrices 
involved are symmetric): 

(2.17) det(/j_fc + fj+k)tf= , det(fj- k - fj+k+2)]^=o, det (/,•_& ± fj+ k +i)]^ =0 - 
They are related to Hankel determinants with symbols on [— 1, 1] by the expressions 

on 2 -2n+2 

(2.18) det(/,_ fc + f J+k )]- k l = —^D%(f(e ie W)/VT^), 
(2-19) det(/,_ fc - f Hk+2 )r k L = ^D»(f(e i9 W)VT^), 

,2 



(2.20) det(/,_ fc + / i+fc+ iffi£o = -^r D n{l{^ {x) )sJY^l 

nn 2 -n /i _ „ 

(2.2i) det(/,-_ fc - Wife = -^r D n(f( el6{x) )]lY^)- 

Asymptotic formulae for Hankel and Toeplitz+Hankel determinants with Fisher-Hartwig 
singularities, whose derivation was based on the above theorems for Toeplitz determinants, 
an asymptotic Riemann-Hilbert-Problem analysis of the polynomials (1.3), and the above 
relations, are presented in [43]. For other asymptotic results, see [16, 17, 8, 13, 36]. 

For related asymptotic results on an important class (see Section 4) of Toeplitz determi- 
nants when the symbol is supported on an arc of the unit circle, see [107, 87, 85, 46]. 

Theorem 2.3 and asymptotic formulae for Hankel and Toeplitz+Hankel determinants find 
applications, e.g., for correlation functions in the XY spin chain in a magnetic field mentioned 
above, in the theory of the impenetrable Bose gas [60, 96], in random matrix conjectures for 
average values of Riemann's ^-function, and Dirichlet's L-functions [82, 66, 35]. 

For a more detailed discussion of the material presented in this section so far, the reader 
is referred to [43]. 

A related area of interest is the asymptotic analysis of Hankel determinants whose symbol 
has Fisher-Hartwig singularities and is supported on the whole real line, or the half-line. 
In particular, in the Gaussian Unitary Ensemble of random matrix theory, the correlation 
function of products of powers of the absolute values of the characteristic polynomial is 
precisely such a Hankel determinant D^(w): namely, the symbol is supported on R and 
given by w(x) = exp(— x 2 ) Y\^ =1 \x — /ij| 2Qj , fij G K, Sfax,- > —1/2. This determinant is 
also related to the 1-dimensional impenetrable Bose gas and conjectures for mean values of 
Riemann's ^-function on the critical line. For a discussion of the results in this area, see 
[86, 70, 71]. For analysis of some other Hankel determinants appearing in random matrix 
models, see [21, 55]. For a recent application of Hankel determinants in the six- vertex model 
see [74, 22, 23, 24, 25]. 

Note that the importance of Fisher-Hartwig singularities appears to stem from the fol- 
lowing feature. The asymptotics of the orthogonal polynomials at the location of such a 
singularity are described by the confluent hypergeometric function [43]. The two indepen- 
dent parameters of such functions are related to the parameters a and /3 of the singularity 
(for /3 — confluent hypergeometric functions reduce to Bessel functions). The location of 
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the singularity corresponds to the single finite branch point of the confluent hypergeomet- 
ric functions. As hypergeometric functions (which depend on 3 parameters) have 2 finite 
branch points, we would not be able to confine ourselves to hypergeometric functions if we 
wanted to consider a singular point which generalizes Fisher-Hartwig in some essential way. 
Roughly speaking, a Fisher-Hartwig singularity is the most general hypergeometric singular 
point. Modifications, of course, are possible: e.g., the end points of the interval [—1, 1] can 
be regarded as modified Fisher-Hartwig singularities for a Hankel determinant with symbol 
1,1] as discussed above. 



on 



3. Transition asymptotics 

A natural question to ask is how the transition between various asymptotic regimes of 
the previous section occurs. Consider once again the 2-spin correlation function for the 2- 
dimensional Ising model discussed above, which is a Toeplitz determinant. As T / T c a 
transition between the Szego asymptotics and the Fisher-Hartwig asymptotics takes place. 
It was first investigated in [112, 93, 101], and the authors found that if T — > T c and n — > oo 
in such a way that x = (T c — T)n is fixed, then the determinant is given in terms of Painleve 
III (reducible to Painleve V) functions. This transition corresponds to the emergence of one 
Fisher-Hartwig singularity with a — 0, /3 — —1/2 at z = 1. The condition that x = (T c —T)n 
is fixed was removed in [38] where uniform asymptotics were obtained for any a, 3la > —1/2, 
(3 G C in terms of Painleve V functions. Namely, consider the following symbol 



(3-1) 



ft(z) 



e t)*+P( z 



e -ty-fi z -^ 



f3 e -in(a+l3) e V(z) 



a ±/3y£ -1,-2, 



where t > is sufficiently small (in the above example of the Ising model, t = const (T c — T)), 
V(z) is analytic in a neighborhood of C, and a, (3 G C with 3?a > — |. The singularities of 
the symbol are at the points e ±l . If t = the symbol possesses a Fisher-Hartwig singularity 
at z = and Theorem 2.2 applies to D n (f ). If t > then ft(z) is analytic in a neighborhood 
of C, and Szego's Theorem 2.1 applies. We have [38] 

Theorem 3.1. Let a, /3 G C with %la > — | and let s$ denote a sector —it/2 + 5 < argx < 
7r/2 — 5, < 5 < 7r/2. Let f t be given by (3.1) and consider the Toeplitz determinants D n (f t ) 
defined by (1.1) corresponding to this symbol. There exists a finite set {x±, . . . ,Xk} G s$ (with 
k = k(a, /3) and Xj = Xj(a, /3) ^ 0) such that there holds the following expansion as n — > oo 
with the error term uniform for < t < to (with to sufficiently small) as long as 2nt remains 
bounded away from the set {x±, . . . , Xk}: 



(3.2) In D n (f t ) = nV + (a + f3)nt + J^ k 



fc=i 



V k -(a + P)- 



-tk 



k 



V- k -(a-P)- 



,-tfe 



k 



. G(l + a + 0)G(l + a-0) „, n , 

G(l + 2a) 



where G(z) is Barnes' G-function, and 
(3.3) 



. 2nt a(x) -a 2 + f3 2 , , 9 n9 , , 

Q(2nt) = / -^ — dx + (a 2 - f3 2 ) ln 2nt. 

x 
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The function o~{x) is a particular solution to the Jimbo-Miwa-Okamoto a-form [76, 77] of 

the Painleve V equation 

(3.4) 

o \ 2 / / \ 2 \ / \ 1 / \ / 

d a\ I da ( da\ da \ I da\ I da \ ( da 



dx 2 ) I dx \dx J dx J \dx ) \dx J \dx 

This solution has the following asymptotics for x > 0: 



a 



2 o2 i a 2 -0 



(5 2 + ^-{x-x 1+2a C(a,(5)}(l + 0(x)), x -+ 0, 2a ^ Z 

(3.5) a(x) = { a 2 -f3 2 + 0(x) + O (x l+2a ) + 0{x l+2a In x), x ->■ 0, 2a £ Z 

r 1+2 " e ~ x n«-mc+n ( l + ° (i)) > x ^ +°°' 

(3.6) o(«,« r(i + « + «r(H»-«r(i-2o) i 



r(l - a + /3)r(l - a - p) r(l + 2a) 2 1 + 2a' 

where T(z) is Euler's V -function. The path of integration in (3.3) is such as to avoid the set 
{xi, . . . , Xk} and is contained within the sector s$. 

Note that (3.4) is the a-form of the Painleve V equation 

/ ~n f l 1 \ o 1 (u-1) 2 ( A B\ Cu u(u + l) 
3.7) u xx =[— + u 2 x - -u x + ^ ^- [Au+-)+ — + D^ -I, 

\2u u-lj x x 2 \ uj x u-1 

with the parameters A,B,C,D given by 

(3.8) A= l -{a-P)\ jB = -I(a + /3) 2 , C = 1 + 20, D = ~. 

The points Xj refer to possible poles of a(x). In the case when a is real and /3 is purely 
imaginary, one can show [38] that a(x) is real analytic for x > 0, and the path of integration 
can therefore be chosen along the real axis. 

If one takes the limit as t — > on the r.h.s. of (3.2), one obtains the correct form of the 
appropriate Fisher-Hartwig asymptotics as given by Theorem 2.2. 

Since the asymptotics of D n (f) are known both at t — and t = t , one obtains an 
amusing identity for the Painleve function a(x): 

G(l + a + l3)G(l + a-P) 

(3.9) fi(+oo) = - In w — ) . 

Methods used in [38] to prove Theorem 3.1 can be adapted to describe other transition 
regimes, e.g., two singularities approaching each other along the unit circle, or emergence 
of an arc on which f(z) = 0. These situations arise for other correlation functions in 
integrable models [61] and appear in the application of random matrix theory to the theory 
of L- functions [35] . 

4. Asymptotics for Fredholm determinants 

We now consider another type of double-scaling limit for Toeplitz determinants which 
yields interesting Fredholm determinants and, after certain analysis, allows us to obtain 
asymptotics of the latter. Note that this approach combined with Riemann-Hilbert-Problem 
techniques allows us to obtain the full asymptotics of these Fredholm determinants including 
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the multiplicative constants which resisted other methods: see [84] for a short review of the 
approach and [85, 41, 42, 9, 46] for details. For a review of other applications of Riemann- 
Hilbert problems to Toeplitz and Fredholm determinants see [48] . For analysis of some other 
Fredholm determinants, see [111, 83, 37] and references in Introduction. 

Let /(e ; n) — 1 on the arc 2s/n < 9 < 2n — 2s/n, < s < n, and f(z; n) = on the rest 
of the unit circle. Then the Fourier coefficients are fo = l — 2s/(nir), fj = — s . , j ^ 0. 
In the limit of growing n and, accordingly, a closing arc, 



vj 



(4.1) ]imD n (f(z;n))=det(I-K smrJ . 



(«) 
where K^ e is the trace-class operator on L 2 (—s, s) with kernel 



sinfa; — v) 
4.2 K aine (x,y) = —f ^. 

In the Gaussian Unitary Ensemble of random matrix theory (and many other random matrix 
ensembles [45, 56]), the Fredholm sine-kernel determinant det(J — K^ e ) describes, in the 
bulk scaling limit, the probability that an interval of length 2s contains no eigenvalues. Of 
interest are the asymptotics of det(J — i*Q|„ e ) when s is large. 

Theorem 4.1. Let K s s ine be the operator acting on L 2 (— s, s), s > 0, with kernel (4-2). Then 
as s — >■ +oo ; 

(4.3) det(J - K { s t) = Csines-'* exp L-*£\ [l + O^ 1 )] , c sme = 2 1 / 12 e^'^\ 

and C'{x) is the derivative of Riemann's zeta function. 

We note that s{d/ds) lndet(7— K^l e ) satisfies [75, 40] a form of the Painleve V equation. In 
particular, this fact enables one to reconstruct the full asymptotic series of the logarithmic 
derivative in the inverse powers of s from the first few terms provided the existence of 
such asymptotic expansion is established. However, the multiplicative constant c S j ne is not 
determined this way. 

Theorem 4.1 was conjectured by Dyson [51] who used, in particular, (4.1) and an earlier 
result of Widom on Toeplitz determinants with a symbol which vanishes on a fixed arc of the 
unit circle [107]. The leading asymptotic term was proved by Widom [106], and the lower- 
order terms apart from c S j ne , in other words the expansion for the derivative (d/ds) lndet(J — 

i*Qi„ e ), by Deift, Its, and Zhou [40] using Riemann-Hilbert methods. Application of a more 
detailed Riemann-Hilbert analysis to Toeplitz determinants allowed the authors in [85, 41] 
to extend the result of Widom [107] to varying arcs, and the relation (4.1) then produced 
the asymptotics of det(J — K^ e ) including c sine , which completed the proof of Theorem 4.1. 
An alternative proof of the theorem was given independently by Ehrhardt [52] who used 
(different) methods of operator theory. 

Recall that f(z;n) was defined above on the arc whose end-points converge to z = 1 as 
n — > oo. We now modify the definition of f(z; n) by placing a Fisher-Hartwig singularity 
at z — 1. Namely, consider the symbol F(z;n) — \z — l\ 2a z^e~ t7Tl3 , z = e te , on the arc 
2s /n < 9 < 2-k — 2s/n, < s < n, and F(z; n) = on the rest of the unit circle. We then 
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obtain 


[46] 










(4.4) 








Um D n (F(z;n)) 
n^oo D n (F(z; oo)) V 




where K^' 


0, S ) . 

IS 


the trace-class operator on L 2 (—s,s) 


with kernel 


(4.5) 




K (f*6 


'{u,v) = 


1 T(l + a + /3)T(l + a- 


• £) A(«)S(u) - A{v)B{u) 


2jci T(l + 2a) 2 


u — V 


where 


















A(x) 


= ^ /2 (a;)|2x| a e- i:c 0(l + a- 


f /3, 1 + 2a,2ix), 



S(x) = gl /2 (x)\2x\ a e ix (j)(l + a-/3, l + 2a, -2ix), 

{-TTi/3 > Q 

' ', «,^C, Ka>-l/2, a±/3^-l,-2,... 

e p , x < 0. 

Here (f)(a,c,z) is the confluent hypergeometric function (see, e.g., [3]) 

v— «, a(a + 1) • • • (a + n — 1) z n 

4.6 a, c, z = 1 + V ) ' K -— -f — . 

v 7 v y *-*i c(c + 1) • • • (c + n - 1) n\ 

The kernel K°^ appears in the representation theory of the infinite-dimensional unitary 
group [29, 26], and the logarithmic derivative (d/ds) lndet(J— fQ^ ) is related to a solution 
of the Painleve V equation. If we set a = (3 = 0, the kernel reduces to the sine-kernel (4.2). 

Theorem 4.2. Let K^ s be the operator acting on L 2 (—s,s), s > 0, with kernel (4-5). 
Then as s —¥ +oo 7 

where G(x) is Barnes' G-function. 

This theorem was proved in [46] using the relation (4.4) and a Riemann-Hilbert analysis. 
The theorem reduces to Theorem 4.1 if a — (3 — (recall that 21nG(l/2) = (1/12) ln2 — 

A particular case of the determinant D n (F(z;n)) with /3 = is related, via a Hankel 
determinant and the formula (2.16), to the following Bessel-kernel determinant det(J — 

^■Besiei) on (0) s )' where the kernel 

( 4 -8) K Bessei( x ' v) = 2(x - y) ' 

and J a (x) is Bessel function. In the Jacobi Unitary Ensemble of random matrix theory (and 
many other ensembles with a so-called hard edge), the Bessel-kernel determinant det(J — 
K-Bessei) describes, in the (left) edge scaling limit, the probability that the interval (0, s) 
contains no eigenvalues. In other words, it describes the distribution of the extreme (smallest) 
eigenvalue. 
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Theorem 4.3. Let K^ s ssel be the operator acting on L 2 (0, s), s > 0, with kernel (4-8) where 

3?a > — 1. Then as s — >■ +oo 7 

(4.9) 

det(/-*t s s L) = c Bessel (a)s- a2 ^ exp (-£ + a^) [l + O^)] , CBeMd ( a ) = ^0- 

These asymptotics were conjectured by Tracy and Widom [102] and proved in [46]. An 
alternative proof of the particular case |5Ra| < 1 is given in [54] by methods of operator 
theory. 

Finally, we consider the case of the so-called Airy kernel. Let w(x) = e~ 4xn be supported 
on J = [0, 1 + s(2n)~ 2 / 3 ] for a fixed s E R, and let D^(w) be the corresponding Hankel 
determinant. Then 

< 4 - io > &, d " + (2^75) = « (' - <» 

where K^- ry is the trace-class operator on L 2 (s, +00) with kernel 
(4.11) tftix,y)= U M M 'to-WM*) . 

Here Ai (x) is the Airy function (see, e.g., [3]). 

In the Gaussian Unitary Ensemble of random matrix theory (and many other ensembles 

with a so-called soft edge), the Airy-kernel determinant det(J — K^ ) describes, in the 
(right) edge scaling limit, the probability that the interval (s, +00) contains no eigenvalues. 
In other words, it describes the distribution of the extreme (largest) eigenvalue. 

Theorem 4.4. Let K^j be the operator acting on L 2 (s, +00), s E R ; with kernel (4. 11). 

Then as s — > — oo ; 

(4.12) 

F TW (s) = det(I-K% ry ) = C ^| S |-V8exp U-^j [l + 0(| S |" 3 / 2 )] , c Mry = 2 1 ^e^- l \ 

The distribution Ftw{s) is known as the Tracy- Widom distribution. Tracy and Widom 
showed that [103] 

(4.13) Ftw( s ) = exp < — / {x — s)u 2 (x)dx 

where u(x) is the Hastings-McLeod solution of the Painleve II equation 

(4.14) u"(x) = xu(x) + 2u 3 (x) , 
specified by the following asymptotic condition: 

(4.15) u(x) ~ Ai (x) as x — > +00. 

The asymptotics of the logarithmic derivative (d /ds) In Fj<\y(s) follow, up to a constant 
(which is in fact zero), from (4.15) and the known asymptotics of the Hastings-McLeod 
solution at —00. The constant c a i ry (as well as CBessei above) was conjectured by Tracy and 
Widom using numerical computations and an analogy with the Dyson formula (4.3). The 
full proof of Theorem 4.4 was given in [42] using (4.10). 
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The determinant det(J — ^llry) a ^ so describes the distribution of the longest increasing 
subsequence of random permutations. Namely, let n = %x% 2 ■ ■ -%n be a permutation in the 
group Sjss of permutations of 1, 2, . . . , N. Then a subsequence i^, ik 2 , ■ ■ ■ ik r , k\ < k 2 < • ■ ■ < 
k r , of 7r is called an increasing subsequence of length r if i^ < ik 2 < • • • < ik r - Let £n(^) 
denote the length of a longest increasing subsequence of tt and let Sn have the uniform 
probability distribution. Then £n{^) is a random variable, and 

(4.16) F TW (s) = lim Prob {tt 6 S N : (l N (n) - 2v / iV)A^ 1/6 < s} 

N— >OD 

This result was obtained by Baik, Deift, and Johansson [6] from the double-scaling limit 
N — y oo, n < N ~ A, of the Toeplitz determinant A nA = D n (exp{\/\(z + z~ v )}). As shown 
earlier by Gessel [63], this determinant is precisely the following generating function: 

oo ^ N 

(4.17) A n A = / u n(N)——T, u n (N) = # (permutations n in Sn with £n(^) < n). 

< * TV' 

7V=0 

An alternative proof of Theorem 4.4 based on this determinant was given in [9]. 

By the Robinson-Schensted-Knuth correspondence (see, e.g., [2]) a permutation tt is related 
to a pair of Young tableaux (of the same shape) of integer plane partitions of N. The number 
In(tt) is the length of the first row of the related tableaux. 

We also note that random permutations are related to last passage percolation and random 
vicious walks [59, 4, 5]. 

There are many related results and extensions of the above results on random partitions 
and permutations, which, in particular, involve asymptotic analysis of special Toeplitz, Han- 
kel, and Toeplitz+Hankel determinants. This large and growing research area has many 
connections to geometry, group representation theory, and integrable models. For details 
and a selection of results, see [6, 78, 79, 28, 95, 7, 8, 73, 57] and references therein. 
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